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ABSTRACT 



If the hot. X-ray emitting gas in rich clusters forms a fair sample of the 
universe (as in Cold Dark Matter (CDM) models), and the universe is at the 
critical density. Ox = 1, then the data appears to imply a baryon fraction, ilb^x 
{ftb,x = derived from X-ray cluster data), larger than that predicted by Big 
Bang Nucleosynthesis (BBN). While various other systematic effects such as 
clumping can lower ftb,x7 in this paper we use an elementary analysis to show 
that a simple admixture of Hot Dark Matter (HDM, low mass neutrinos) with 
CDM to yield mixed dark matter shifts ilb^x down so that significant overlap 
with Ofc from BBN can occur for Ho^75 km/sec/Mpc, even without invoking 
the possible aforementioned effects. The overlap interval is slightly larger for 
lower mass neutrinos since fewer cluster on the scale of the hot X-ray gas. We 
illustrate this result quantitatively in terms of a simple isothermal model. More 
realistic velocity dispersion profiles, with less centrally-peaked density profiles, 
imply that fewer neutrinos are trapped and, thus, further increase the interval 
of overlap. However, we also note that if future observations of light element 
abundances find that Ofc/i^^O.OlS, the range of concordance in this simple mixed 
dark matter model vanishes. 
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1. Introduction 

Big Bang Nucleosynthesis (BBN) has repeatedly been found to be in agreement (within 
the tolerances of the measurements) with observed light element abundances. Along with 
these observations, BBN predicts a baryon density, 

0. 007/^-2^06^0. (1) 

(Copi, Schramm, and Turner 1995 a, b, c) where h = i7o/(100 km/sec/Mpc). Recent work 
(Barlett et al. 1994, White et al. 1993) using the ROSAT and ASCA X-ray sateUites 
has demonstrated that rich clusters of galaxies contain significant amounts of hot X-ray 
emitting gas. If one assumes that the gas is gravitationally confined, then the temperature 
of the gas can be related to the gravitational potential and, hence, the total mass. My, of 
the X-ray emitting region. The intensity of the emission is related to the density of the gas 
(baryons) in the emitting regions. The mass of the galaxies within the cluster is computed 
based upon photometric luminosities of the galaxies (e.g., Godwin and Peach 1977) and an 
empirical mass-luminosity relation (van der Marel 1991). If these clusters are fair samples 
of the universe, then the ratio of the hot gas mass plus the luminous mass in the galaxies, 
Mfc, to the total mass, Mt is a measure of ^h,x^ the value of derived from X-ray cluster 
data, [Mb/Mr = 06,^/Ot). 

White et al. (1993) show that the current data (using the simplest "P* order" analysis) 
imply 

0.005 + 0.03/^-3/2^^^0.013 + 0.07/^-3/2^ ^2) 

where first term refers to the galaxies and the second term refers to the X-ray gas. Fig. 1 
(adapted from Gott et al. 1974, Copi and Schramm 1995) shows the comparison of 
from BBN (Eqn. 1) and ilb^x from M^/Mt (Eqn. 2) assuming that the clustered gas is 
representative of the baryon density in an Ox = 1 universe. Notice that there is only overlap 
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for Ho^bb km/sec/Mpc. Of course, if Ox < 1 or if the clusters are not a fair sample, then 
there is no conflict for larger values of Hq. Furthermore, systematics (c.f. Bird, Mushotzky, 
and Metzler 1995) seem to be capable of shifting M^/Mt primarily to lower, not higher, 
values. For example, density inhomogeneities in the hot gas would produce an anomolously 
high implied gas mass. Also, if significant magnetic fields were present in the cluster such 
that the magnetic pressure was non-negligible compared to the thermal pressure, then 
Mt would be larger than the value derived by considering the thermal pressure alone. 
Furthermore, weak lensing estimates of Mt (Squires et al. 1995) seem to yield slightly 
larger values than those from the X-ray data, which could further lower M^/Mx. In this 
paper, however, we will accept the M^/Mt taken from the simplest "P* order" analysis of 
the X-ray data but will argue that a Hot Dark Matter (HDM) component (neutrinos) can 
allow for a greater overlap. Obviously in a pure HDM model with structure formed by the 
wakes of cosmic strings (Brandenburger 1994, Stebbins 1995) clusters are unlikely to be 
fair samples. However, in this paper we will show that even in standard, infiation-inspired 
models an HDM admixture can increase the interval of overlap. 

In pure Cold Dark Matter (CDM) models, the scale of clusters (i?^Mpc) have obtained 
a fair sample of the universe and, hence, 

Ofc,, = nrMb/Mr. (3) 

Since the most interesting case is Ox = 1, we will focus our attention on this option. When 
one switches to a mixed dark matter model, with HDM as well as CDM, Ox in Fqn. 3 
becomes the cosmological density of all of the matter which clumps on the scale of the 
cluster, iOduster = bilciusten whcrc ftduster = 0.2 it 0.1 (Fabcr and Gallagher 1979) and 
b ~ few is a bias factor which need not be the same as the galaxy formation bias factor 
used in some CDM models of structure formation. Thus, 

/^duster = Ofc -|- OcDM + / ' ^HDM, (4) 
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where / is the fraction of HDM which clumps within the cluster. Mixed dark matter 
(HDM + CDM) models (Occhinero 1985) have been shown to be a viable way to enable 
COBE-normalized CDM-like models to fit galaxy correlation data, for ftcDM = 0.8 with 
flj = 1 (Primack et al. 1995, Davis 1995, Ma and Bertschinger 1994). 



2. Low Mass Neutrinos 



Neutrinos are the best candidate for the HDM. The sum of the masses of all of the 
species of neutrinos (and anti-neutrinos) yields (Lee and Weinberg 1977, Cowsik and 
McClelland 1972, Marx and Szalay 1972) 

^^^^ = 9#^- 

The fraction of neutrinos which will contribute to the cluster is / = £hdm I ^(hhCDM)^ where 
e^; is the ratio of the density of component x within the cluster to the cosmological density 
of the same component. In the spirit of the Thomas-Fermi corrections used in models of 
degenerate electrons within a Coulomb potential (Salpeter 1961), we can write £hdm for a 
single neutrino mass (which could be degenerate and shared by multiple neutrino species) 
locally as a ratio of Fermi integrals, £hdm = -^(?7ij/</)(r)/A;rj,)//(0), where (/)(r) = —GM[r)/r^ 
Tj, = 1.9 K, and 

is the relativistic Fermi distribution function. (A relativistic distribution is employed since, 
due to the low cross section for neutrino interactions, the distribution has not changed 
since the neutrino temperature exceeded about 1 MeV (Lee and Weinberg 1977, Tremaine 
and Gunn 1979).) We can generalize to neutrino species with different masses by replacing 
f-^HDM in Fqn. 4 with fj-^HDM{j), where fj = f{m^.) and ilnDMU) = 2m^J{92h^ eV), 
with the factor of 2 arising from summing over Uj and Uj. For the remainder of this paper 
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we will consider only the case of a single, possibly degenerate, neutrino mass, as this will 
serve to suitably illustrate the features of this model. 

We defined e such that Xduster = £(bkCDM){^b + ^cdm) + £hdm^hdm, where 



where XComa ~ 280 (e.g.. White et al. 1993). If < £hdm ^ Xduster as one would expect, 
then < / < 1. Fig. 2 shows how / depends on the strength of the cluster potential, (t^^), 
with the cluster radius fixed at 1.5h~^ Mpc, the standard Abell radius. Fig. 3 gives the 
fraction of trapped neutrinos as a function of the radius of the cluster with the potential 
strength fixed at (1800 km s~^)^, the value adopted by White et al. (1993). The cluster is 
parameterized in terms of these two quantities {{v'^) and R). Notice from Fqn. 6 that £hdm 
is a function of (u^) = (/> oc M/ R but not R alone and that Xduster is a function of (u^) and 
R] Xduster M / R^ OC {v"^) j R^ . Thus, Fig. 2 illustrates the effects of changing £hdm and 
Xduster siucc (u^) is an independent variable in both parameters. However, in Fig. 3 only 
Xduster varics since £hdm is not a function of R alone. Fig. 4 gives / as a function of niu. 
Note that / increases with niu since more massive neutrinos are more likely to clump on 
the scale of rich clusters. 

In computing £hdm we adopt an isothermal velocity dispersion; therefore, £hdm is 
constant throughout the cluster. This assumption is quite good for our purposes, since even 
fairly large variations in the actual mass distribution of the cluster have little effect upon 
the derived fraction of trapped neutrinos. See Fig. 2. Furthermore, an isothermal model 
is the most conservative choice in that it over-estimates the number of neutrinos bound 
to the cluster. In an isothermal model, p oc . Any model in which the total density 
profile is less centrally peaked (Bird, Mushotzky, and Metzler 1995, Girardi et al. 1995) 
will have |</>(r)| < |</>(-R)| everywhere within the cluster. Thus, / will be smaller than in the 
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isothermal case. (See Eqns. 6 & 7.) 

Fig. 5 shows the overlap in as derived from BBN (Eqn. 1) and our mixed dark 
matter interpretation of the Coma cluster data (Eqns. 2, 3, & 4) for ftcDM = 0.8. It is clear 
that this overlap increases with decreasing neutrino mass, since fewer neutrinos will clump 
on the scale of the cluster. However, Ox = 1, ^b,x is fairly insensitive to for a particular 
value of ilcDM- To see this, we can expand ilb^x in powers of /; from Eqns. 3 & 4 we have 



for ftcDM = 0.8 and h = 0.5. Eor this value of /i, Eqn. 5 specifies a maximum neutrino mass 
(one uu-paii) of ~ 2 eV. Eig. 4 shows that /^0.005 for this range of neutrino masses. 
Thus, changing the neutrino mass, or similarly the number of massive neutrino species, will 
result in a change in ilb^x on the order of several hundredths of a percent. 

Although we have shown that the X-ray cluster data can be made to be consistent 
with the predictions of Big Bang Nucleosynthesis for Ho^75 km/s/Mpc, we note that this 
concordance depends sensitively upon the upper limit of O^/i^ derived from the observed 
light element abundances. If future measurements of light element abundances find that 
Ofc/i^ lies near the present upper limit ijlj^h^ 0.024), then the consistency of X-ray 
cluster data and BBN will be retained. However, if O^/i^ is found to be less than the value 
representing the intersection of the lower O^^^; curve and the line h = 0.38 ijlj^h^ 0.018), 
there will no longer be a range of concordance. 

The maximum value of Hq = 100 /i km/s/Mpc consistent with BBN and these cluster 
data is given by 



Eig. 6 illustrates the upper limit of the range of consistency for /i as a function of ^cdm 
for these observations for a range of ^cdm and for various numbers of massive neutrino 



Ofc,^ oc 1 



Ofc -|- OcDM 



/ + ...- 1 + 0.2/, 
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species. Since oc VIhdm^ the limiting value of h becomes insensitive to the multiplicity 
of the mass degeneracy for larger values of VIcdm^ where ftnoM tends toward zero. This 
figure also illustrates this decreasing importance of the actual number of massive neutrino 
species with increasing ilcDM- 

3. Conclusions 

Structure formation in a Cold + Hot Dark Matter model is initiated by perturbations 
in the CDM field. The neutrinos are trapped in the primordial fiuctuations and thereby 
accrete onto the new structure. Our analytical model illustrates that neutrinos contribute 
to the cluster halo. Although the model is simplistic, its results are believable. We find 
that neutrinos are only a small component of the mass of the cluster. Since the neutrinos 
are such a small component of the mass of the cluster, Ox is effectively represented by the 
(smaller) coduster as described in the text. Thus, we have shown that in the framework of a 
cold + hot dark matter model with massive neutrinos these data are consistent with BBN, 
Ox = 1, and Ho^75 km/sec/Mpc. In such a mixed dark matter model the more "radical" 
options of Fukgita, Hogan, and Peebles (1993) for resolving this issue become unnecessary. 
As long as observed light element abundances allow Ofc/i^^O.OlS, there will be a range of 
concordance between X-ray cluster data and BBN. This critical value of il^h^ will increase 
if low values of Hq are ruled out. However, as mentioned earlier, numerous systematic 
effects could shift the cluster data, lowering this critical value and increasing the range of 
concordance toward larger values of Hq. 

We acknowledge useful discussions with Mike Turner, Craig Copi, and Evalyn Gates. 
While this work was in the final stages of preparation we received a preprint from Kofman 
et al. (1995) in which quantitatively similar results were obtained. However, their paper 
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de-emphasized the increased overlap due to their focus on higher values of Hq. This work 
was supported in part by NSF, NASA, and DOE at The University of Chicago and by 
NASA at Fermilab. 
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Fig. 1. — Values of allowed by Big Bang Nucleosynthesis and implied by X-ray cluster 
observations for a CDM model with Ox = 1- Concordance between the two methods of 
determining requires that i7o^55 km/sec/Mpc. Also shown are bounds based upon the 
age of globular clusters (to > 10 Gyr) and Type la SN determinations of Hq [Hq > 38). 

Fig. 2. — The fraction of trapped neutrinos as a function of virial velocity dispersion 
(normalized to the value for the Coma cluster). We have assumed ftcDM = 0.8 and = 0.04 
(consistent with BBN for Hq 50 km/sec/Mpc). Thus the sum of all neutrino and anti- 
neutrino masses is approximately ISh^ eV. 

Fig. 3. — The fraction of trapped neutrinos as a function of cluster radius (normalized to the 
radius for the Coma cluster) for a fixed (u^) = {v'^)coma- See Fig. 2 for details concerning 
OcDM, Ofc, and m^. 
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Fig. 4. — The fraction of neutrinos of a given mass (in eV) which are trapped within the 
Coma cluster for an ftcDM = 0.8, 0^ = 0.04 model. Note that an ftnoM ~ 0.2 universe will 
not allow a uu-paii with individual masses in excess of ~ 5.5 eV for h^O.S. Thus, m^k^b.b eV 
requires that Oc/^m^O.S. 



Fig. 5. — Values of concordant with Big Bang Nucleosynthesis and X-ray cluster 
observations. The cluster model assumes that ftcDM = 0.8, Ox = 1, /i = 0.5, and = 2 eV, 
consistent with a single massive uu-paii. Values of consistent with both methods require 
Ho^75 km/s/Mpc. See Fig. 1 for details about the other features of this plot. 



Fig. 6. — Maximum values of h from concordance of Big Bang Nucleosynthesis and "P* 
order" X-ray cluster observations as a function of ilcDM- Note that these curves do not 
represent particular values of since the total mass in neutrinos is a function of h and 
^HDM = 1 — (Ofc -|- (Icdm)- The inset shows the dependence of the concordance interval 
upon the number of massive neutrino species. The three curves converge for larger values 
of ftcDM since oc ftnoM tends toward zero. Thus, for reasonable values of ilcDM^ the 
number of massive neutrino species is unimportant, provided that there is at least one. 



